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We describe searches for B meson decays to the charmless vector-vector final states ojK* and top 
in 89 million BB pairs produced in e + e~ annihilation at y/s = 10.58 GeV. We measure the following 
branching fractions in units of 10" 6 : B(B° luK*°) = 3.4±^;| ± 0.4 (< 6.0), B(B + uK* + ) = 
3.5l^ ± 0.7 (< 7.4), B(B° -> up ) = 0.6^ ± 0.4 (< 3.3), and B(B+ -* top + ) = 12.6±^ ± 1.6. 
The first error quoted is statistical, the second systematic, and the upper limits are defined at 
90% confidence level. For B + — > Lop + we also measure the longitudinal spin alignment fraction 
f L = 0.88lg:i5 ± °- 03 and charge asymmetry A ch = (5 ± 26 ± 2)%. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The BABAR 0,SI3 and Belle experiments have re- 
ported observations of B meson decays to most of the 
charge states of two-body combinations of p and K* 
mesons. Here we present the results of searches for decays 
to final states with an uj meson plus a K* or p meson. 

The decays B — > V4V2, where V\ and V2 are spin-one 
mesons, proceed through a combination of S-, P-, and D- 
wave amplitudes, or in the helicity basis by amplitudes 
A , A±i. The subscripts give the helicities of the vector 
mesons. The longitudinal Aq amplitude is a linear combi- 
nation of S- and D-wave (CP even), while the transverse 
A±i contain all three partial waves 

The spins of the vector mesons are analyzed by their 
decays into pseudoscalar mesons. We define the helicity 
frame for each vector meson as its rest frame with po- 
lar axis along the direction of the boost from the B rest 
frame. Because of the limited statistics of the present 
study we integrate over the azimuthal angles, exploiting 
the uniform azimuthal acceptance. The angular distri- 
bution is 



T d cos Videos 62 

I - h) sin 2 9 1 sin 2 6 2 + f L cos 2 B x cos 2 2 J , 

where 9k is the decay angle (defined below) in the V% 
helicity frame, and f L = |A>l7(l^+i | 2 + |A>| 2 + \A^\ 2 ) 
is the fraction of the longitudinal spin component. Mea- 
surements in B — ► pp decays find values close to fr, = 1 
IH 0- Q > consistent with expectation [(| , while smaller 
values found in B — ► cf)K * Q are theoretically puzzling. 
Since the different angular distributions for the longi- 
tudinal and transverse components lead to different ac- 
ceptances, the /^-dependence needs to be considered in 
searches as well as in detailed studies of B — > V1V2 de- 
cays. 

In the decays of B + (or the flavor-definite B° — > 
K*°X, K*° —>■ K + ir~) and their charge conjugates it is 
also of interest to measure the direct- CP- violating charge 
asymmetry A ch = (F~ - r+)/(r~ + T+) in the rates 
r± = T{B ± for each final state /±. 

In this study we search for the four decays B° — > uiK*°, 
B+ -> luK*+, B° ojp°, and B+ -> ojp+ @. A previ- 
ously published search by CLEO || established 90% con- 




00 ( f ) 

FIG. 1: Representative Feynman diagrams for B — > top and 
B — > ujK* decays: (a) external tree; (b, c) CKM-favored glu- 
onic penguins; (d) CKM-suppressed penguin; (e, f) destruc- 
tively interfering color-suppressed trees. 



fidence level (C.L.) upper limits on the branching frac- 
tions of (23, 87, 11, and 61)xl0 -6 , respectively Because 
these charmless B decays involve couplings with small 
CKM mixing matrix elements, several amplitudes poten- 
tially contribute with similar strengths, as indicated in 
Fig. ^ The B + modes receive contributions from ex- 
ternal tree, color-suppressed tree, and gluonic penguin 
amplitudes, with the external tree (a) favored for B + — > 
ujp + , and the penguin (b) strongly favored by CKM cou- 
plings for B + — > u)K* + . For the B° modes there are no 
external tree contributions, and again, for B° — ► ujK* 
the penguin (c) is CKM-favored. For B° — ► ujp° the 
color-suppressed tree amplitudes (e, f) almost cancel jjjj 
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TABLE I: Selection requirements on the invariant mass (in 
MeV) and decay angle of B-daughter resonances. 



State 



mv. mass 



decay angle 



p° 510 < m(7T7r) < 1060 -0.85 < cos6> < 0.85 

p+ 470 < m(7T7r) < 1070 -0.6 < cos 9 < 0.85 

K"°+^,K*+ Qv+ 755 < m(Kir) < 1035 -0.85 < cos<9 < 1.0 



K 



K° s 



755 < m{Kix) < 1035 
735 < m(7T7r7r) < 825 
120 < m(ry) < iso 
488 < m(7T7r) < 508 



-0.6 < cos6> < 1.0 
< I cos (91 < 1.0 



because of the different isospins of the final-state mesons, 
leaving only a Cabibbo-suppressed penguin (d). Weak 
exchange and annihilation amplitudes are expected to be 
negligible. 

Theoretical estimates of the branching fractions for 
vector-vector decays include those based on isospin re- 
lations among various modes E3 > effective Hamiltonians 
with factorization and specific i?-to-light-mcson form fac- 
tors GH El El El, and QCD factorization The es- 
timated branching fractions lie in the range < 10 -6 (for 
B" -> top ) to 20 x 1(T 6 (for B+ -> ujp+). 

The results presented here are based on data col- 
lected with the BABAR detector at the PEP-II asym- 
metric e + e~ collider E3 located at the Stanford Lin- 
ear Accelerator Center. An integrated luminosity of 
81.9 fb _1 , corresponding to 88.9 ± 1.0 million BB pairs, 
was recorded at the T(4S) resonance (center-of-mass en- 
ergy y/s = 10.58 GeV). 

Charged particles from the e + e~ interactions are de- 
tected, and their momenta measured, by a combination 
of five layers of double-sided silicon microstrip detectors 
surrounded by a 40-layer drift chamber, both operating 
in the 1.5-T magnetic field of a superconducting solenoid. 
We identify photons and electrons using a CsI(Tl) elec- 
tromagnetic calorimeter (EMC). Further charged parti- 
cle identification (PID) is provided by the average energy 
loss (dE/dx) in the tracking devices and by an internally 
reflecting ring imaging Cherenkov detector (DIRC) cov- 
ering the central region. 

We reconstruct the i?-daughter candidates through 



, p+ -> 7r+vr u , K*» 
K+^{K*+ + ^ ), K*+ 



their decays p° — ► n 
K+n-(K*^ +w _), K*+ 

K°n + (Kt± ,), us -> tt+tt-tt , tt° -> 77, and K° 

7r + 7T _ . Table [i] lists the requirements on the invariant 
mass of these particles' final states. For the p, K* , and us 
invariant masses these requirements are set loose enough 
to include sidebands, as these mass values are treated 
as observables in the maximum-likelihood fit described 
below. For candidates we further require the three- 
dimensional flight distance from the event primary vertex 
to be greater than three times its uncertainty. Secondary 
pions and kaons in p, K* , and us candidates are rejected 



if their DIRC, dE/dx, and EMC PID signature satisfies 
tight consistency with protons or electrons, and the kaons 
(pions) must (must not) have a kaon signature. 

Table H| also gives the restrictions on the K* and p he- 
licity angles 9 made to avoid regions of rapid acceptance 
variation or combinatorial background from soft parti- 
cles. We define 9 as the angle relative to the helicity axis 
of: the normal to the decay plane for us, the positively- 
charged (or only charged) daughter momentum for p, and 
the daughter kaon momentum for K*. 

A -B-meson candidate is characterized kinemat- 
ically by the energy-substituted mass Toes = 
[{^ s + Po ' Pb) 2 /Eq — p%]* and energy difference 
AE = Eg — \\fs, where the subscripts and B refer to 
the initial T(AS) and to the B candidate, respectively, 
and the asterisk denotes the T(4S) frame. The reso- 
lution on AE (m ES ) is about 30 MeV (3.0 MeV). We 
require \AE\ < 0.2 GeV and 5.20 < m ES < 5.29 GeV. 
The average number of candidates found per selected 
event is in the range 1.15 to 1.2, depending on the final 
state. We choose the candidate with the smallest value 
of a x 2 constructed from the deviations of the daughter 
resonance masses from their expected values. 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e + e~ — > qq events (q = 
u,d,s,c). We reduce these by selecting on the angle #t 
between the thrust axis of the B candidate in the Y(4S) 
frame and that of the rest of the charged tracks and neu- 
tral calorimeter clusters in the event. The distribution 
of |cos6*t| is sharply peaked near 1.0 for combinations 
drawn from jet-like qq pairs, and nearly uniform for B- 
meson decays. The requirements, which optimize the ex- 
pected signal yield relative to its background-dominated 
statistical error, are |cos6>t| < 0.8 for the K* modes 
and |cos#t| < 0.65 for the p modes. In the maximum- 
likelihood fit we also use a Fisher discriminant T E3 that 
combines four variables defined in the T(4S") frame: the 
angles with respect to the beam axis of the B momen- 
tum and B thrust axis, and the zeroth and second angu- 
lar moments L0.2 of the energy flow about the B thrust 
axis. The moments are defined by Lj — J2iPi x l c °s#i| J , 
where 0, is the angle with respect to the B thrust axis 
of track or neutral cluster i, pt is its momentum, and the 
sum excludes the B candidate daughters. 

From Monte Carlo (MC) simulation E3 we estimate 
the residual charmless BB background to be 0.1% or less 
of the total sample in all cases. To allow for contributions 
possibly missing in the simulation we include a compo- 
nent for these in the fit described below, with a yield free 
to vary. 

We obtain yields, fa, and A c h from extended un- 
binned maximum-likelihood fits with input observables 
AE, toes, J 7 ) and for vector meson k the mass m/j 
and helicity- frame decay angle 9k- For each event i and 
hypothesis j (signal, continuum background, BB back- 
ground) we define the probability density function (PDF) 
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TABLE II: Signal yield Y and bias Yo with their statistical uncertainties, detection efficiency e, daughter branching fraction 
product, significance S (with systematic uncertainties included), measured branching fraction, and 90% C.L. upper limit for 
each mode. The number of produced B mesons is (88.9 ± 1.0) x 10 6 . 



Mode 


Y 

(events) 


Yo 
(events) 


(%) 


(%) 


S 
(*) 


B 

(io- 6 ) 


B U.L. 

(io- 6 ) 


u>K*° 




3.2 ±1.1 


13.2 


59 


2.2 


3.4+^g ± 0.4 


6.0 


uK k+*» 

u>K*+ 


11-63:5 

tr 4+6.O 


2.9 ±1.1 
-0.1 ±0.8 


13.3 
6.7 


20 
30 


1.3 
1.4 
1.9 


3.9+3^ ±0.9 
3.1+J-4 ±0.9 
3.5+1* ± 0.7 


7.4 




4 3+11.0 


-0.5 ± 1.0 


10.5 


89 


0.4 


0.6+^1 ± 0.4 


3.3 




rrj .-.+ 18.5 

°'-'-16.5 


-1.3 ±2.8 


5.4 


89 


4.7 


12. 6+I 3 ± 1.6 





V) = V 3 {m^)V 3 {^E l )V 3 {r)V 3 {m^ 2 ,e\^ 2 ). (2) 

We check for correlations in the background observables 
beyond those contained in this PDF and find them to be 
small. For the signal component, we correct for the ef- 
fect of neglected correlations (see below). The likelihood 
function is 



e -(E^) Jt„ 



(3) 



where Yj is the yield of events of hypothesis j found by 
maximizing C, and N is the number of events in the 
sample. 

The PDF factor for the resonances in the signal takes 
the form Vi. sis (m\)V 2 , s i g {m 2 )Q{0\,e 2 ) with Q given by 
Eq^modified to account for detector acceptance. For qq 
background it is given for each resonance independently 
by V q7} (ml,dl) = V pk (mi)V pk (ei) +V c (mi)V c (0i), dis- 
tinguishing between true resonance (Ppk) and combi- 
natorial components (T^c)- For the BB background 
we take all four mass and helicity angle observables to 
be independent. The other PDF forms are: sum of 
two Gaussians for "PsigC^Es), 'Psig(AE), and the peak- 
ing components of Vj(m^); a conjunction of two gaus- 
sians with different widths below and above the peak 
for Vj{T)\ and linear or quadratic dependences for AE, 
rrik , and helicity cosines for qq combinatorial background. 
The qq background in hies is described by the function 
x\J 1 — x 2 exp [— £(1 — x 2 )] , with x = 2niEs/V^ an d pa- 
rameter £. 

For the signal and BB background components we de- 
termine the PDF parameters from simulation. We study 
large control samples of B decays to charmed final states 
of similar topology to verify the simulated resolutions in 
AE and mES> adjusting the PDFs to account for any 
differences found. For the continuum background we use 
(m,ES,AE) sideband data to obtain initial values, before 
applying the fit to data in the signal region, and ulti- 
mately leave them free to vary in the final fit. 



Free parameters of the fit include signal and back- 
ground yields, background PDF parameters, and for the 
mode for which we find a significant signal, and the 
signal and background charge asymmetries. For the fits 
without significant signal we vary /l by hand, basing the 
upper limit on the fit with = 0.9, a choice that is con- 
sistent with a priori expectations and gives a conservative 
efficiency estimate. The free background-PDF parame- 
ters are £ for mES, slope for AE, area and slope of the 
combinatorial component for m&, and the peak position 
and lower and upper width parameters for T. 

We evaluate possible biases from our neglect of corre- 
lations among discriminating variables in the PDFs by 
fitting ensembles of simulated experiments into which we 
have embedded the expected number of signal events ran- 
domly extracted from the fully simulated MC samples. 
We give in Table [H] the values found for bias for each 
mode. Events from a weighted mixture of simulated BB 
background decays are included, and so the bias we mea- 
sure includes the effect of crossfeed from these modes. 

In Tabled we show for each decay mode the measured 
branching fraction together with the quantities entering 
into its computation and with its uncertainty and signifi- 
cance. The statistical error on the signal yield or branch- 
ing fraction, fi,, and A c h is taken as the change in the 
central value when the quantity —2 ln£ increases by one 
unit from its minimum value. The significance is taken 
as the square root of the difference between the value of 
— 21n£ (with systematic uncertainties included) for zero 
signal and the value at its minimum. For all modes ex- 
cept B + — ► up + we quote a 90% C.L. upper limit, taken 
to be the branching fraction below which lies 90% of the 
total of the likelihood integral in the positive branch- 
ing fraction region. In calculating branching fractions 
we assume that the decay rates of the T(4S) to B + B~ 
and B°B a are equal. For decays with K* + , we combine 
results from the two K* decay channels by adding the 
values of — 21n£, taking into account the correlated and 
uncorrelated systematic errors. 

We present in Fig. |2 the data and PDFs projected 
onto mES and AE, for subsamples enriched with a mode- 
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FIG. 2: Projections of jties (left) and AE (right) with a cut 
on the per-event signal/total likelihood ratio for (a, b) B° — > 
loK*°, (c, d) B+ -» loK*+, (e, f) B° -> lo P °, and (g, h) B+ -> 
ujp + . The solid (dashed) curve gives the total (background) 
PDF, computed without the variable plotted, and projected 
onto the same subspace as the data. 
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FIG. 3: Background-subtracted projections of helicity cosines 
for (a) lo and (b) p + from the fit for B + — > Lop + . 



dependent threshold requirement on the ratio of signal to 
total likelihood (computed without the PDF associated 
with the variable plotted) chosen to optimise the signif- 
icance of signal in the resulting subsample. Fig. [3] gives 
background-subtracted projections onto the helicity an- 
gle cosines for B + — > uip + corresponding to the fit result 
/l = 0.88^0 15 i 0.03; the dominance of the term pro- 
portional to jh in Eq. ^is evident. 

The branching fraction value B given in Table ^ for 
B + — > u>p + comes from a direct fit with the free pa- 
rameters B and /l, as well as A c h- This choice exploits 



the feature that B is less correlated with than either 
the yield or efficiency taken separately. The behavior of 
-2\nC(f L ,B) is shown in Fig.H 
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FIG. 4: Distribution of x 2 = -21n£(/ L ,6) for B+ lu P + . 
The solid dot gives the central value; curves give the contours 
in 1-sigma steps (A^/x^ = 1) out to four sigma. 



Most of the systematic uncertainties on the branching 
fractions arising from lack of knowledge of the PDFs have 
been included in the statistical error since most back- 
ground parameters are free in the fit. For the signal, 
the uncertainties in PDF parameters are estimated from 
the consistency of fits to MC and data in control modes. 
Varying the signal-PDF parameters within these errors, 
we estimate yield uncertainties of 1-4 events, depending 
on the mode. The uncertainty in the fit-bias correction is 
taken to be half of the correction itself. Similarly we esti- 
mate the uncertainty from modeling the BB backgrounds 
by taking half of the contribution of that component to 
the fitted signal yield. These additive systematic errors 
are dominant for the modes with little or no signal yield. 
We have also considered backgrounds from B decays to 
the same ultimate final state as the signal. The helic- 
ity angle restrctions given in Table [I] suppress ujtt or u)K 
subsystems in the region of known resonances. For the 
B° — > cop and B° — > uiK*° upper limits, inclusion of 
the helicity-angle PDF with fixed Jl = 0.9 reduces to 
a negligible level the effect of interferences with possible 
S-wa.ve 7t — (it, K) states. 

Uncertainties in our knowledge of the efficiency, found 
from auxiliary studies, include 0.8% x N t , 2.5% x 7V 7 , 
and 4% for a K® decay, where Nt and 7V 7 are the num- 
ber of tracks and photons, respectively, in the B candi- 
date. Our estimate of the i3-production systematic error 
is 1.1%. Published data |2(j provide the uncertainties in 
the .B-daughter product branching fractions (1%). The 
uncertainties in the efficiency from the event selection are 
1-3% for the requirement on cos^t and 1% for PID for 
the modes with a charged kaon. The dependence of effi- 
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ciency on fa causes uncertainties of 2-6% in the B — > loK 
and B — > up measurements. 

The 0.03 systematic error on fjj for B + — > a>p + comes 
from imperfect representation of correlations in the PDF, 
and is estimated from fits to fully simulated MC sam- 
ples. From several large inclusive kaon and i?-dccay sam- 
ples, we find a systematic uncertainty for A c h of 2% due 
mainly to the dependence of reconstruction efficiency on 
the charge of the p-daughter charged pion. The value of 

A q c q h = ± °- 7 ) % that we nnd for the background in 

the B + — > ujp + fit provides confirmation of this estimate. 

In summary, we have performed searches for the pre- 
viously undetected decays B° — > wi\T*°, 5+ — > 



<° and B+ 



We observe B + 



ujp 



with 



a significance of 4.7 <7, and establish improved 90% C.L. 
upper limits for the other modes: 



(3.4±i:| ±0.4 (< 6.0)) x 10~ 6 
(3.5±|jj±0.7 (< 7.4)) x 10~ 6 
(0.6±J:f ± 0.4 (< 3.3)) x 10~ 6 
(12-6±1;| ± 1.6) x 10~ 6 , 

where the first error quoted is statistical and the second 
systematic. For B + — > uip + we also measure the longitu- 
dinal polarization fraction 



B(B° -> ljK*°) 
B{B+ -> wiT+) 

S(£ + -► ^p+) 



0.88^5 ±0.03 



and charge asymmetry 



A ch = (5 ±26 ±2)%. 

We find that the longitudinal spin alignment is domi- 
nant, as for the pp modes yj, |2j, |3j, |4j . The central value 
of the branching fraction for B + — > is about half 
of those found for £? + — * p+p and B° — * p + p~ . All of 
our branching fraction results are in general agreement 
within errors with the theoretical estimates. 
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